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AN EXPLICIT VAN DER CORPUT ESTIMATE FOR £(1/2 +it) 


GHAITH A. HIARY 


Abstract. An explicit estimate for the Riemann zeta function on the critical 
line is derived using the van der Corput method. An explicit van der Corput 
lemma is presented. 


1. Introduction 

The Riemann zeta function is defined for s = cr + it by 

oo 

( 1 ) C(s):=^n" s , cr > 1. 

n= 1 

It can be analytically continued everywhere except for a simple pole at s = 1. A 
well-known problem in number theory is to bound the growth rate of zeta on the 
critical line a = 1/2. This problem has led to deep ideas in the theory of exponential 
sums. In particular, the method of exponent pairs (see [ 6 ], [14]; page 116]), and the 
Bombieri-Iwaniec method mm- 

Since |£(cr + *t)| = |£(cr — it)\, we can suppose that t > 0. Starting with the series 
Q, it follows by the Euler-Maclaurin formula that £(1/2 + it) <C t 1 / 2 for large 
enough t. This can be improved substantially by appealing to the Riemann-Siegel 
formula, which gives £(1/2 + it) <C t 1 / 4 . The Weyl-Hardy-Littlewood method (see 
m section 5.3]) detects a certain amount of cancellation in the main sum of the 
Riemann-Siegel formula, further improving the bound to £(l/2 + *t) <C t 1 / 6 log 3 / 2 t. 
The van der Corput method (e.g. [Mill]) removes an extra y/\ogt factor, which 
sharpens the estimate to 

( 2 ) £( 1/2 + it) <C £ 1 / 6 log£. 

The exponent 1/6 in estimate (0 is hard to improve. The sharpest result so far 
is due to Huxley [9], who proved that £(1/2 + it) <C t 32 / 205 log 7 1 for some constant 

7 . Note 32/205 = 0.15609_ (See also the recent result in [3].) Assuming the 

Riemann hypothesis, one can show that £(1/2 + it) -C exp(Alogf/log log t) for 
some constant A; see [14] §14.4], Hence, £(1/2 + it) is conjectured to grow slower 
than any fixed power of t. 

In this article, we obtain an explicit bound of the van der Corput type 0. That 
is, we compute constants C\ and C 2 such that |£(l/2 + if)| < C/t 1 / 6 logt for t > C 2 . 
This is part of our work in progress about subconvexity bounds for zeta. We are 
also motivated by the following computational considerations where the simplicity 
of an explicit van der Corput estimate is particularly attractive. 
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Specifically, in numerical tests of the growth rate of £(1/2 + it), it is necessary 
to have explicit bounds. For one usually cannot distinguish a small power of t from 
a logarithm or a subexponential factor unless t is prohibitively large. An explicit 
bound offers an unconditional measure against which one can compare large values 
of |C(l/2 + it )| found by special numerical searches, such as those in [7]. 

Another motivation comes from the algorithms in [8] which employ a Taylor ex¬ 
pansion to express zeta as a sum of low degree exponential sums and their deriva¬ 
tives. The size of the remainder term in these algorithms is bounded by the highest 
order derivative used and the maximal size of a certain subdivision of the main sum 
of zeta. In this context, an improved and simple method to obtain an explicit bound 
is of value as it enables reducing the number of derivatives needed to guarantee a 
given error tolerance. This in turn will improve the running time appreciably. 

Explicit bounds for zeta were obtained by Cheng and Graham |4j, and recently 
by Platt and Trudgian jT2] who proved that 

(3) |£(l/2-Mt)| < 0.732f 1/6 logf, (t > 2). 

We improve the leading constant in <(3|) by 14%. 

Theorem 1.1. If 0 < t < 3, then |£(l/2 + it)\ < 1.461. If t> 3, then 

(4) |£(l/2 + it)\ < 0.63t 1 ^ 6 logf. 


A numerical computation reveals that max t >3 |£(l/2 + it)\/t 1 ^ 6 logt > 0.507. 
So, a priori, the leading constant in (j4|) cannot be reduced below 0.507 without an 
additional assumption on the size of t. In fact, if one employs the Riemann-Siegel 
Lehman bound |£(l/2 + it)\ < 4(f/(27T) 1 / 4 — 2.08 from Lemma 12731 for any range of 
t at all, then the leading constant cannot break 0.541, as a numerical computation 
shows that 


( 5 ) 


mm ■ 

t> 3 


4(f/(2 7 r) 1 / 4 - 2.08 


> 0.541. 


t 1 / 6 log t 

Our proof of Theorem 11.11 uses a subdivision of the Riemann-Siegel main sum 
different than H ng, giving rise to a simpler optimization problem. In particular, 
we divide the main sum into short pieces of length ss t 1 / 3 (which appears to be a 
new, and natural, subdivision; c.f. Weyl’s method in in Section 5.3] which divides 
the main sum into pieces of length < f 1 / 6 ). This subdivision enables better control 
of the oscillations in each piece via Lemma [L2l since the range of f"'{x) for us will 
be more restricted. 


Lemma 1.2. Let /( x) be a real-valued, function with three continuous derivatives 
on [AT +1, N + L\. Suppose there are W > 1 and A > 1 such that < |/ ,,, (*)| < ^ 
for N+l<x<N + L. Ifr]>0, then 


( 6 ) 


N+L 

E 

n=N -\-1 


a 2irif{n) 


< (LW~ 1/3 + rj)(aL + /3W 2/3 ), 


where 


( 7 ) 


a := a(W,\,ri) = ± + ^Vv + W~^ + ^ 

8 ■= 8(W n) = — — H_ 

P ■ P{ ,V) 15y/rj + Wi/ 3 ' 


A 

1R2/3’ 
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Remarks. Lemma Q is proved in 0 This lemma is an explicit version of the 
process AB in the method of exponent pairs. The condition W > 1 in the lemma 
can be relaxed to W > 2/n. The number 77 will be chosen of size ss 1 in our 
application, Theorem 11.11 

The proof of Theorem 11.11 completely overlooks cancellation among the ~ f 1 / 6 
pieces where Lemma 11.21 is applied. This is the main source of inefficiency in our 
proof. As far as we know, there is no definitive mechanism to take advantage of this 
cancellation. Instead, one works with longer pieces, and tries to prove cancellation 
within each one. 

Of course, the bound 0 is asymptotically far from the truth. And even for 
moderately large values of t, this bound is still probably a substantial overestimate. 
Evidence for this comes from computations by J. W. Bober and the author, some of 
which are summarized in ]7j- In these computations, several hundred large values of 
zeta were recorded by computing £(1/2 + it) at certain special points. The largest 
value found this way was 

(8) |C(l/2 + LZo) | = 16244.86526..., 

where T 0 = 39246764589894309155251169284104.050622 ..., so T 0 « 3.9 x 10 31 . (To 
our knowledge, © is the largest value of zeta computed so far.) In comparison, 
the bound 0 gives |C(l/2 + iTo)\ < 8448744, which is about 521 times larger 
than 0. Alternatively, one can use the van der Corput bound (I2U1) directly. This 
gives a bound that is about 507 times larger than our computed value. It should be 
stressed, though, that @ was found by searching a thin set of special points, not by 
an exhaustive search, and so only provides a lower bound for max tg [ 0 i T 0 ] IC(l/2+*t)|- 

2. Proof of Theorem 11.11 
For the remainder of the paper, we set 

(9) c 0 := 0.63, t 0 ■= 9.3 x 10 7 . 

Proof. We consider t over three ranges. In the range 0 < t < 200, we use the 
computational bound from Lemma 12.21 In the intermediate range 200 < t < to we 
use the Riemann-Siegel-Lehman bound supplied by Lemma 12.31 And in the last 
range t > to we use the van der Corput bound supplied by Lemma 12.41 □ 

To handle the intermediate and last ranges of t, we rely on Lemma 12.11 which is 
consequence of the Riemann-Siegel formula. This lemma requires t > 200, which 
is the reason we treat the initial range t £ [ 0 , 200 ] separately. 

Lemma 2.1. If t > 200 and n\ = [\/t/(27r)J, then 

Wl 

(10) |£(l/2 + it) | < 2| Y, n~ 1/2+it \ + U(t). 

n= 1 

where K(t) := lASt" 1 / 4 + 0.127t“ 3 / 4 . 

Proof. The lemma gives a small improvement on [12] Lemma 3]. The improvement 
is that the inequality for 7 Z(t) is tighter if t > 200. The lemma is stated separately 
for emphasis, as it is an essential first step in all that follows. 
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We apply the triangle inequality to the Riemann-Siegel formula in [5] page 9] to 
obtain 


( 11 ) IC(l/2 + it)\ < 2 


E 

n= 1 


cos(tlogn — 9(t)) 


y/n 


( 2tt 


1/4 


+ P ( “ ) + l^o(t)|, 


where n± = |_\/t/ 27 rj, |i?o(OI < 0.127t 3 / 4 for t > 200 , and 

n(^ 2 

( 12 ) 


B = max 

0<z<l 


cos | (z 2 + 3/4) 


cos 7 TZ 


By® page 65], we have (3 < 0.93 (more precisely, /? = cos(7t/8)). The lemma now 
follows from the inequality [3(2i r ) 1 / 4 < 1.48, and using 


(13) 


m 

v 

cos(t logn — 9(t)) 

1 

< 

ni At log n—i6(t) 

1 

711 it log 7i 

n— 1 

yjn 

- 2 

^ Vn 

71=1 

+ 2 

71=1 


E’ 

71—1 


— 1/2 +it 


□ 


In choosing the cutoff point for the intermediate range, we computed an approx¬ 
imation, call it fi, for the largest t > 200 such that the Riemann-Siegel- Lehman 
bound (flHl) beats the van der Corput bound (l20l) . This gave t\ ~ to- To decide 
the leading constant in Theorem ll.il we computed an approximation, call it ci, for 
the smallest c > 0 such that the Riemann-Siegel-Lehman bound evaluated at t\ is 
smaller than ct\' 6 \ogt\. This gave c\ « cq. 

The computational bound in Lemma |2.2I was verified using interval arithmetic, 
which is also the method used in [ 121 . 

Lemma 2.2 (Computational bound). If 0 < t < 3, then |£(l/2 + it)\ < 1.461. If 
3 < t < 200 , then |C(l/2 + it)\ < cq t 1 / 6 logt. 

Proof. We implemented the Eulcr-Maclaurin formula (see (TTJ [13]) using interval 
arithmetic. We remark that a substantial loss is incurred in the upper bound (and 
computational speed) due to the use of interval arithmetic, but this loss is still 
tolerable for our purposes. 

We used a main sum of length [2 1] terms in the Euler-Maclaurin formula with 
a single correction term. Given an interval [ao, 6 o]> our program computed an 
enveloping interval I out = [oi, &i] such that /2 + it)\ : t £ [a, b]} C I ou t- Hence, 

(14) max |£(l/2 + it)\ < max |t| = b\. 

t^[ao,bo] t£lout 

With this in mind, we partitioned the interval [3, 200] into consecutive subintervals 
Iq := [q/Q i (q + 1)/Q], where q = 3 Q, ..., 200Q — 1 and Q = 2 7 . For each I qi our 
program returned an enveloping interval I q , 0 ut which we used in (1141) to verify that 

maxtgj, |C(l/2 + it)\ ^ 

( ’ (g/Q)V 6 i 0 g(g/Q) - C °- 

The estimate m held for all relevant q , in fact, with a smaller constant of 0.595. 

To verify the bound |C(l/2 + it)\ < 1.461 for t £ [0,3], one could prove that 
|C(l/2 + it)\, 0 < t < 3, attains its maximum at t = 0. However, this seemed 
unduly complicated. Instead, we fell back on our interval arithmetic program after 
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modifying some its parameters. Specifically, we used the Euler-Maclaurin formula 
with a main sum of |~60(t + 1)] terms, which is much longer than before. We kept 
a single correction term, and also used a finer partition with Q = 2 14 . The longer 
main sum for 0 < t < 3 ensured that the error in the the Euler-Maclaurin formula 
was sufficiently small. Given this, we were able to verify the claimed bound. □ 

Lemma 2.3 (Riemann-Siegel Lehman bound). Ift> 200, then 

4t 4 / 4 

(16) |C(l/2 + it)|<—-^-2.08. 

In particular, for 200 < t < t 0 , we have |£(l/2 + it)\ < Cot 1 / 6 \ogt. 

Proof. This lemma gives a small improvement on the Lehman bound in US Lemma 
2]. The improvement is in the term —2.08, which is significant if t is not too large. 

We start with the bound furnished by the Riemann-Siegel Lemma [Til To bound 
the main sum there, we employ the estimate 



= 2 —= + 4 y/nf — 4\/5. 

, V n 

n= 1 


Note that n\ = | \\ftj (27r)J > 5 for t > 200 , so the integration step in (fTTl) makes 
sense. Also, we have 2 X]n=i n ~ 1 ^ 2 ~ 4\/5 < —2.48. To bound the remainder term 
7 Z(t) in the Riemann-Siegel Lemma 12.II we employ the estimate 

(18) K{t) < 1.48(200)- 1/4 + 0.127(200)- 3/4 < 0.4, (f > 200). 

The first part of the lemma now follows on substituting CD and m back into 
m, and noting that y/nf < (t/ 27r) 1 ^ 4 . 

To prove the second part of the lemma, we use Mathematica to verify that there 
is no solution to the equation 

4+ 1 / 4 

(19) WP* - 2- 08 = coi 1/6l °gL (200 <t< t 0 ), 

and that the l.h.s of the equation is smaller than the r.h.s. at t = 200 , and therefore 
throughout the range t £ [ 200 , to]- □ 

Lemma 2.4 (Van der Corput bound). If t > to, then 

( 20 ) |C(l/2 + *t)| < ait 1 / 6 logt + a 2 t 1 / 6 + a 3 

where cq := 0.6058490462530, a 2 := 0.5743984045897, and a 3 := -2.884626766806. 
In particular, fort > to, we have |£(l/2 + it)\ < cot 1 / 6 logt. 

Proof. We plan to divide the main sum in the Riemann-Siegel Lemma 12.11 into 
pieces of length « t 1 / 3 , then apply the van der Corput Lemma fL2l to each piece. 
To this end, let K = [t 1 / 3 ] and R = [_ n\/K \, where, as before, m = \_\Jtf 27tJ . 
Here, K is the length of each piece (except possibly the last one, which can be 
shorter), and R+ 1 is the total number of pieces. 
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The remainder term TZ(t) := 1.48t 1 / 4 + 0.127t 3 / 4 in the Riemann-Siegel 
Lemma m satisfies 1Z < TZ(to) for t>to- Thus, carrying out the aforementioned 
subdivision, and using the triangle inequality, we obtain 


( 21 ) 


r 0 K — l 


R-l (r+l)if—1 


IC(l/2 


+ ,:i )is2 E th +2 S| E 

n=1 v r=rn n=rt 




n ^ git log n 


RK 


sjn 


—To n=rK 

K{to). 


e it log n 



Here, we trivially estimated the part of the main sum with n < r^K, where ro is a 
positive integer to be chosen later. Also, we used R > 0, which is due to f > to- 
To bound the first sum in pTj) , we note that r$K > \rot 0 ' ]. So, proceeding as 

in urn we obtain 


( 22 ) 

where 

(23) 


r 0 K -1 

2 E 


n =1 



< 4 \Jr$K +T(r 0 ,t 0 ), 


RotJ /3 l-l - 

Ar 0 , t 0 ) := 2 —j=- 4 V rr-oto /3 l - 1- 

i V ^ 

r).— l v 


To bound the remaining sums in (EH), we use partial summation [14] (5.2.1)]. Put 
together, if we let 


(24) 


S:= 2 




max 


V r K A<if 

r=rg 


A —1 

^ ( git log(rA'+fc) 
fe =0 


then we obtain 


(25) |C(l/2 + tt)| <5 + 4 x /ro^ + I(r 0 ,t 0 )+^(to). 

In order to bound the inner sum of S, we employ the van der Corput Lemma ll.2l 
Setting 


(26) 
then 

(27) f"\x) = 
So, on defining 

(28) W:=W r = 


t 


f(x) ■= — log (rK + x), 


ir(rK + a;) 3 ’ 
7 T (r + 1) 3 AT 3 


and noting that A < A', we obtain 
(29) 


W s s W’ 


(0<x<A-l), 

(0 < x < A - 1). 

x ._ x __ ( r + l ) 3 

A . A r i 

(0 <y< K). 


We apply Lemma fT72l with L = K, a r := a{W r , X r ,rj), (3 r := /3(W r ,rj), and with 
r] > 0 to be chosen later. This yields 

R I - 

(30) S < 2 Y —== Ja r K 2 W^ 1/3 + r]a r K + /3 r I<Wr /3 + ril3 r W? /3 . 

_ v rK 
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We factor out K 2 W r ^ 3 = n ifs[r+i) from under the square-root. This gives 

(31) 
where 

„ , mrW r 1/3 , (3 r W? /3 , V p r W r 

(32) B r := a r H-—-1-—-h 

K K 

Since K < f 1 / 3 + 1 and R < t 1 / 6 /y/V/K, we have 

(33) Wr <T:{R + l)\l+t-^f<^j=, 

In addition, K > t 1 / 3 , so put together we obtain 

rja r p f5 r p 2 


1/3 


2^/6 « 
* - 7T 1 / 6 ^ 


/ B r 
r(r + 1) ’ 


K 2 


p := 


(34) 


J3 r ^ Oi r + 


( 1 + i? 

VPrP 3 


1 + 


**)• 


V / 2'7T 1 / 6 i 1 / 6 27T 1 / 3 2 3 / 2 v / 7rf 1 / 6 

At this point, we make several observations. First, lF r is monotonically increasing 
in r, and A r is monotonically decreasing in r. Hence, by the definitions of a r and 
Pr, we see that they monotonically decrease with r. So a r < a ro and /3 r < /3 ro for 
r > ro. Second, we fix r] = rj o := (75/64) 2 / 3 , which is to help balance the leading 
two terms in the formula for a r in 0. Third, if we define 

\Ao/(27t) - 1 


(35) 


Ro := 


C +1 


- 1 


then R> Ro, and therefore 

(38) P<Po:= (l + ^)(l + ^a)- 

Assembling these estimates into (13T1) . and using 

R 


(37) 


V'' — - < log ■ 

™ y/r(r + 1 ) r 0 - 1 


R 


together with the inequality t > to to bound the denominators in (13411 from below, 
we obtain 


(38) 


S < 


2 i 1 / 6 


VO&r 0 pO 


Pr 0 Po 


~ ttVB 


VoPraPo log _ R 


r ° v^TrVe^/ 6 2 tt 1 / 3 2 3 / 2 v ^4 /6 "'° r 'o- 1 ' 

We substitute (1551) back into (1551) . use the inequality W ro > 7r(ro + l) 3 to simplify 
the bounds for a ro and fi ro , and choose ro = 5 which is suggested by numerical 
experimentation. Then, we combine the resulting expression with the bounds 


(39) 4 < 4y/r 0 (l + t~ 1/3 )t 1/6 , R < t 1/6 /V2 n, 

and numerically evaluate the resulting constants using Mathematica. On comple¬ 
tion, this yields the first part of the lemma. 

To prove the second part of the lemma, denote the r.h.s. of CU1) by (*), and 
consider the equation (*) = cot^logt. Using Mathematica, we find that there is 
no solution t > to to this equation, and that (*) is smaller than the Coi 1//6 logt at 
t = to, and therefore for all t > to- □ 
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Remarks. The reason we restrict t > to in Lemma 12.41 is that the Riemann-Siegel 
Lehman bound from Lemma 12.31 is tighter for t < to- So we may as well take t > to, 
which gives marginally better constants. 


3. Proof of Lemma [T~2l 


Proof. We use the Weyl-van der Corput Lemma in |4j Lemma 5], but in the more 
precise form presented at the bottom of page 1273. Also, we incorporate a refine¬ 
ment pointed out by Platt and Trudgian in m Lemma 2] that allows writing the 
leading term in (l40ll as L + M— 1 instead of L + M. Put together, if M is a positive 
integer, then 


(40) 
where 

(41) 


N+L 

E 


g 2 iri/(n) 


n=N -\-1 


2 

<(L + M 



sI(>-b)iwi). 


N-\-L—m 


S'm(L) '■= J2 


g 27ri(/(r+m)-/(r)) 


r=N +1 


Henceforth, we may assume that m < L and L > 1. Otherwise, the sum S' m is 
empty, and so it does not contribute to the upper bound in flUl) . 

Let g(x) := /(x + m) — f(x), where N+l<x<N+L— m. Then, g"(x) = 
/"( x + m) — f"(x). Therefore, by the mean-value theorem, g"{x) = mf '"(£,*) 
for some f* £ (x, x + m). Since we assumed that 1 < m < L and L > 1 then 
(x,x + m) C [N + 1, N + L], So, given our bound on we deduce that 

(42) ^ < \g"{x)\ < (N + l <x< N + L-m). 

Applying the van der Corput Lemma in [4j Lemma 3] to S' m (L) thus yields 


(43) 


1 0 r ml / 8A Lyjm/W , 3A Lm , 8 sJWjm , 0 
\S m (L)\ < - + -^r- + 5 + 3. 


We substitute (TS1) into 6 DD , then execute the summation over m using the esti¬ 
mates that appear after j4] Lemma 7]; namely, 


(44) 




4M 3 / 2 

15 


E( 


1 - 


m \ 

MJ 



4 VM 
3 


as well as the Euler-Maclaurin summation estimates 

M T, r M 


(45) 


£(‘-s) 


< 


M 


m =1 




m \ 


MJ 


m < 


M 2 

IT 


From this, we conclude that 


(46) 


M 


< 


32A LM 3 / 2 
75 VW 


XLM 2 i 32 VWM | 3 M 
2 W + 15 + ~2~ 


In particular, substituting fi6l) back into we obtain 


N+L 

E 


g 2ir if(n) 


n=N+l 


2 

< {L + M- 1) 



64AL [M XLM 
75 V W + W 


64 

15 
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We choose M = [ 77 W 1 / 3 ] for some free parameter 77 > 0 that can be optimized 
(usually, 77 will be around 1). This choice is in order to balance the first two terms 
on the r.h.s. as they typically dominate in our application. So, now, appealing to 
the inequality rjW 1 ^ < M < rjW 1 ^ 3 + 1, we deduce that 


(47) 


N+L 

E 


e 2ir if(n) 


n=N-\-l 




64A L I r/W 1 / 3 + 1 
75 V W 


XLinW 1 / 3 + 1) 64 / W 

+ w + i5\j^w^ + 3 


We factor out LW 1 ^ 3 from the first three terms in the second bracket, and W 1 ^ 3 
from the last two terms. This gives 


N+L 

E 

n=N-\-l 


p 2n if(n) 


<(L + ? 7 W 1 / 3 )-Erf- + ^07 + W~ 1 / 3 


+ 


W 1 / 3 V 77 75 

\{r] + W- 1 / 3 ) 


)+{L + r,W 1/3 )W 1/3 + 3r 1/s ) . 


W V 3 


150? 


Last, the lemma follows on recalling the definitions of a and /3 in 0 - 


□ 
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